Grain boundary melting in ice 
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We describe an optical scattering study of grain boundary premelting in water ice. 
Ubiquitous long ranged attractive polarization forces act to suppress grain boundary 
melting whereas repulsive forces originating in screened Coulomb interactions and 
classical colligative effects enhance it. The liquid enhancing effects can be manipu- 
lated by adding dopant ions to the system. For all measured grain boundaries this 
leads to increasing premelted film thickness with increasing electrolyte concentration. 
Although we understand that the interfacial surface charge densities q s and solute 
concentrations Ci can potentially dominate the film thickness, we can not directly 
measure them. Therefore, as a framework for interpreting the data we consider two 
appropriate q s dependent limits; one is dominated by the colligative effect and one is 
dominated by electrostatic interactions. 
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I. INTRODUCTION 



Surface and interfacial premelting have been studied for the principal facets of ice and 
for ice crystals in contact with a wide range of other materials 1 . Despite the ubiquity of 
polycrystalline materials in laboratory and natural environmental settings, and the basic 
relevance of grain boundaries in effective medium properties, direct measurements of grain 
boundary (GB) premelting for systems in thermodynamic equilibrium have not been made. 
Crystallography insures that the mismatch between two grains is characterized by molecular 
scale disorder, the structure of which, as the melting temperature is approached from below, 
is basic to the edifice of premelting. However, in general liquid-like grain boundary structure 
is a controversial topic with a complex relationship between film thickness, temperature 
and chemical composition 2 ' 3 . For example, differing interpretations arise from studies in 
metals 1 '', colloidal crystals 6 , and molecular solids such as benzene'. Simulations suggest 
complete 8 and partial melting 51 , or the formation of a third orientation of the solid phase 10 . 
A similar diversity of behavior is observed between grains in multicomponent systems 2,11 ' 12 . 



The study of ice is compelling for many reasons 1 , a few of which we mention here. 
First, it exhibits the same class of phase transitions found in more simply bound matter. 
Second, it can be held near the pure bulk melting point T m with relative ease. Third, 
the results are of immediate relevance to geophysical phenomena. In order to probe the 
grain boundary between two ice crystals a light scattering apparatus was constructed to 
nucleate ice bicrystals, control their growth and to expose them to varying levels of ions 
using dissolved salt 13 . As thermodynamic parameters are varied, laser light is reflected 
from a single grain boundary and the intensity of the reflected signal is interpreted using a 
theory that incorporates the optical anisotropy of the bounding crystals 11 . In that theory 
we analyzed the reflection and transmission of plane waves by an isotropic layer (water) 
sandwiched between two uniaxial (ice) crystals of arbitrary orientation. The experimental 
geometry is depicted schematically in Fig. 1 and here we briefly summarize the approach 
detailed in Thomson, Wettlaufer, and Wilen 13 . 



2 




FIG. 1. Schematic of the optical bench geometry as described in 13 . A laser is spatially filtered and 
frequency chopped before the polarization is set and the beam is focused onto the grain boundary 
(GB). The reflected signal is refocused onto a calibrated photo detector and data is output to a 
PC at selected time intervals. 

II. EXPERIMENTAL METHOD 

Our thermally insulated ice growth apparatus is used to nucleate and grow ice in contact 
with a solution of monovalent electrolyte impurity (NaCl). The apparatus is mounted onto 
rotation and translation stages and placed in the path of a laser beam on an optical bench. 
The incident beam is spatially filtered, frequency chopped, and polarized before it is focused 
onto the grain boundary and the reflected signal collected and analyzed. In order to limit 
absorption and scattering within the ice, use of a low power 2.3 mW, 632.8 nm He-Ne laser 
insured that less than 0.02% of the beam was attenuated within the ice itself. Hence, any 
heating effect was insignificant and below the resolution of our thermometry. The reflected 
signal is refocused onto a photovoltaic detector whose output goes to a lock-in amplifier 
(SR830), and frequency-locks the measured signal to the incident beam. The frequency- 
locked signal and temperature data are read by a computer at specified time intervals. 
Bicrystals, with different crystallographic mismatch, are grown and we continuously collected 
the reflected intensity data as a function of the temperature and impurity concentration. 
From the measured voltage at the lock-in an intensity ratio is calculated after applying 
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a series of corrections. The adjustments include accounting for the response of the photo 
detector, the measurement algorithm of the lock-in, and attenuation of the beam intensity by 
the optical elements of the system. Finally, the intensity ratio is converted to a film thickness 
using the theoretical model described in Thomson, Wilen, and Wettlaufer 1 1 . Consistent with 
other measurements 15 ' 1(1 it is assumed that the index of refraction of the film is the same 
as bulk water (n w = 1.33) which has weak thermal 1 ' and solutal ls dependence. To obtain 
a measurement of d from the intensity ratio we average many measurements taken at fixed 
bulk composition and temperature giving a combined maximum error in the measurement 
and intensity conversion of ± 3%. Therefore, the reflected light intensity measurements 
allow us to determine GB thicknesses as a function of temperature, impurity concentration, 
and independently measured crystal orientation. 

The key features of the data can be summarized as follows: The range of experimental 
salt concentrations Cj of the bulk liquid in contact with our ice grain boundaries is 0-5 
psu, and the range of observed grain boundary of thicknesses is 1-8 nm, and we have no 
direct measurement of the surface charge density q s . At zero concentration of electrolyte, a 
nonzero GB thickness is observed, with different values for GBs of various misorientations. 
As the solute (NaCl) concentration is increased, there is a trend of increasing GB thickness. 
In general there is some more pronounced noise in this trend at lower solute concentrations. 



III. THEORY 

A theoretical treatment of the role of impurities in interfacial premelting aides in the inter- 
pretation of the data. The approach predicts that in the presence of soluble electrolyte impu- 
rities GB premelting in ice may be macroscopic very near the bulk melting temperature 20,21 . 
A GB liquid film of thickness d arises from the competition between attractive van der 
Waals interactions, colligative effects and repulsive interactions of electrostatic origin; the 
presence of a surface charge q s is mediated by ions within the liquid layer. Minimizing a free 
energy that combines all of these effects 19 ' 20 ' 22 provides the equilibrium film thickness d at 
an impure grain boundary held at a temperature T as follows; 

(1) 

Here, pi and q m are the liquid density and the latent heat of fusion, AT = T m — T is 
the undercooling relative to the melting temperature of a pure ice-water system, iVj is the 
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FIG. 2. Experimental grain boundary thickness d(nm) as a function of NaCl concentration Cj(psu) 
determined as described in the text. Colors represent different grain boundary orientations deter- 
mined from an independent measurement via polarimetry and a Schmidt plot as described in 13 . 
The orange and red data correspond to the grain boundaries with the largest lattice mismatch, 
and hence the largest magnitude of interfacial disorder at contact in the limit of zero impurity 
concentration. The blue data show a small angle twist grain boundary and hence the weakest 
interfacial disorder at contact. Dark green (olive) points correspond to the same grain boundary 
as the concentration is cycled up (cycled down and back up) showing hysteresis. The lines are 
fits for the small q s colligative limit described by Eq. (4). Finally, the vertical arrow corresponds 
to the upper bound of the surface charge density dependent critical composition above which the 
theory of predicts that d increases with composition. The upper and lower bounds depend on 
q s , the former corresponding to ~ 0.2 C m -2 (the vertical arrow) and the latter to ~ 0.02 C m -2 , 
which would be denoted by a vertical line on the origin. As seen in Fig. 3 when q s changes through 
this range so too does the compositional range of non-monotonic variations in d. The average 
colligatively shifted undercooling for these data is approximately 1.5 K. 



number of moles of electrolyte per unit area, Ah is the Hamaker constant, e and e Q are the 
relative and free space permittivities, k _1 = \J ^z^n- ls ^ ne Debye length, and R g , k B , Na 
and e are the ideal gas and Boltzmann constants, Avogadro's number and the elementary 
charge respectively. Although in general the frequency dependence of the dielectric response 
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of ice/water systems can lead to important retarded potential effects 1 that can in principle 
influence a grain boundary film 20 , such effects are presently too subtle to be probed. More 
importantly, the Hamaker constant of relevance in Eq. (1) is Ah = 3.3 x I0~ 22 , which is 
calculated from full frequency dependent dispersion force theory 19 , renders this attractive- 
film suppressing-term negligible for the range of film thicknesses observed experimentally. 
The question here is what limits of Eq. 1 can explain the experimental data? There are two 
limits of the theory that are relevant to consider. 

We take one limit where q s very small, and there is very slow equilibration of composition 
along the GB. Here, the effect is principally colligative and hence d increases with C\. The 
small q s limit of Eq. (f) and vanishing dispersion forces is represented by Eq. (29) of 
Hansen-Goos and Wettlaufer 19 . As shown in Fig. 3, the theory applies to the case when 
the solute along the grain boundary (or grain boundary network in bulk ice) is isolated 
and the absolute undercooling is specified. The total number of moles, and hence areal 
concentration, iVj, of solute is conserved. This is relevant, and has been applied, to glaciers 
and polar ice where the diffusion of impurities is slow enough that the grain boundaries in 
the interior may be considered to be isolated. 

In the second limit q s plays a larger role and we fit its dependence on Cj in a manner 
consistent with Eq. (I). Here, d also increases with d in the GB, which is equilibrated to 
the bulk fluid reservoir. The dominating effect is the screened Coulomb repulsion between 
the two sides of the GB. We discuss these limits in turn. 

A. Small q s , Slow Equilibration, Locally Colligative Limit 

A motivation for this limit is that data are qualitatively consistent with an assumption 
that the actual concentration Aj along the grain boundary is proportional to the bulk con- 
centration of the reservoir into which the ice grows. Thus, one may write Aj = C^eff where 
£ e ff is a length scale that relates the two quantities. Physically, one can think of this length 
scale as being related to a "funneling" effect from the grain boundary groove as the ice lens 
is grown into the bulk solution, trapping the total solute across this effective length into the 
grain boundary. Hence we are then led to consider the possibility that, on the time scale 
over which the measurement is performed, the bulk concentration in the GB remains greater 
than that of the reservoir either (a) due to slow diffusion, possibly influenced by finite size 
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effects, or (b) due to effects that lower the chemical potential along the grain boundary in 
a way not captured in the current theory, but along the lines discussed in 23 . 

This limit is plotted in Fig. 3 where we see behavior similar to the experimental fluc- 
tuations at low Ni and the experimental trend of increasing film thickness with the areal 
impurity level. There is a sensitive dependence of the film thickness on q s at very low dopant 
levels, where the slope changes sign below a threshold value N t ^ ~ 8 q s /xM m~ 2 . However, 
when Ni > iV t h the value of d always increases with iVj. 

To be more quantitative, if we postulate the above, then we can express the colligative 
limit of Eq. (1) as follows: 
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where AT' is the undercooling with respect to the solutally depressed melting temperature. 
Upon rearranging, we find 
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Finally, we further assume that there is a minimum thickness d Q to the grain boundary due 
to short range effects and hence we can write 

piq m AT' 
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where d Q will depend on the particular grain boundary mismatch. The straight lines in Fig. 2 
are the fits using Eq. (4) for which we find an average £ e g of approximately 30 nm and used 
the average AT' for each experiment. 

B. Large q s , Fast Equilibration, Electrostatic Limit 



An experimental motivation for this limit is that the measured GB always extends to 
the bulk fluid reservoir outside the ice, and the temperature is measured relative to that at 
the ice-solution edge. Thus, we consider the case when the solute in the GB is diffusionally 
connected to an external reservoir of specified bulk solute concentration. In equilibrium, 
then, one might expect the volume concentration of solute in the grain boundary to be equal 
to that of the reservoir. In this limit, one may replace Ni/d by C$ in Eq. (1). Rearranging 
terms yields 
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FIG. 3. The theoretical dependence of the film thickness d on impurity level Ni in moles of 
electrolyte per unit area (here written as mol m -2 ) for the range of surface charge density q s 
from 0.02 to 0.22 C m~ 2 . The undercooling AT is fixed at the average experimental value of 
approximately unity. The surface charge q s increases from 0.02 in the lowest curve to 0.22 in the 
upper solid curve in increments of 0.04 C m~ 2 and the calculation becomes meaningless below a 
cut off of approximately the molecular scale. The essential points are (a) that d depends sensitively 
on q s at very low dopant levels, (b) for TV, > iVth — 8 q s /zmol m -2 , d will always increase with Ni 
(for any undercooling) and (c) over a range of undercooling there is a range of q s at small impurity 
concentrations where d decreases as iVj increases, as discussed in Hansen-Goos and Wettlaufer 19 . 

As already noted, we have no experimental measurements of the interfacial surface charge 
and we have no direct measure of the composition within the GB. Thus, we take as an 
ansatz that q s is a function of the concentration Cj. The experimental data shown in Figs. 
2 and 5 demonstrate a substantial scatter but nonetheless motivate one to consider a linear 
dependence of the grain boundary thickness d on the impurity concentration Cj. Thus, for an 
undercooling of 1.5 K we take d = 3 nm for Cj = 0.5psu and d = 4 nm for Ci = 1.6psu and 
linearly extrapolate. Now, neglecting the temperature dependence of q s over the relatively 
small experimental temperature ranges, we use Eq. (5) to calculate q s (Ci), the result of 
which is shown in Fig. 4, and we use this in Eq. (5) to determine ci as a function of Ci for 
undercoolings from 0.5 to 2.0 K. We compare these theoretical results with the experimental 
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FIG. 4. The surface charge density q s as a function of the concentration C, of the bulk fluid 
reservoir assumed to be in equilibrium with the GB. 

data in Fig. 5, in which we see that beyond Cj ~ lpsu the general linear increase in the GB 
thickness with concentration is captured. The inter-experimental variability at low values of 
Ci is again seen here; depending on the undercooling d can either increase or decrease with 
Ci. 

IV. DISCUSSION AND CONCLUSION 

We do not have a testable microscopic description relating the measured crystallographic 
orientations of the ice crystals denning the GB to the nature of the disorder at short range. 
Such a theory may facilitate some interpretation of the variation of the observed intercepts 
d in Fig. 2, but we note that they range from about 4 to 24 molecular layers, and bulk fluid 
properties have been demonstrated in subnanometer scale water films 15 ' 16 . As seen in Figs. 
3 and 5 for a wide range of values of the surface charge density, the regions with JVj < N t h 
exhibit the largest variation in film thickness with q s . Hence, in this context we interpret 
the observed fluctuations in d ~ d Q as being associated with the combined sensitivity of the 
magnitude of d on q s and of the sign and magnitude of the slope of the film thickness with 
impurity level. 

Each GB is grown with different thermal and solutal histories, although the intergranular 
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FIG. 5. Theoretical prediction of grain boundary thickness d(nm) as a function of Cj(psu) deter- 
mined as described in the text surrounding Eq. (5) for the large q s , electrostatic limit. From the 
lower line to the upper line the undercoolings are 2.0, 1.5, 1.0 and 0.5 K. The data are replotted 
with the same symbols as in Fig. 2 but here to avoid confusion we leave off the fits from Eq. (4). 
We note that the data are all ostensibly at the same undercooling of about 1.5K. It is not our 
intent to fit the data from any given GB in this figure, but simply to show the overall behavior of 
this limit of theory and the range over which d increases with d. 

ions can exchange and relax with the bulk reservoir, the fluctuations associated with these 
histories leave unique ionic concentrations behind, bound to differently charged surfaces. As 
described above, it is evident that the effects associated with additional sensitivities when 
N{ < iV t h can amplify such history dependence at low dopant levels. Because we do not have 
intra-GB film measurements of ion concentration, to speculate further than is warranted by 
direct evidence risks interpretation that cannot be tested quantitatively with our method. 

These measurements provide essential information for the effective medium properties of 
ice polycrystals. For example, the interfaces between grains provide a ready pathway for the 
transport of soluble impurities and isotopes, which in the case of Earth's ice sheets act as high 
resolution proxies for past climates. Indeed, it has been shown that the liquid veins where 
three grains abut and the nodes where four or more grains terminate can act as the principal 
conduits for transport of such proxies 23 . However, we know that the total intergrain surface 
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area is dominated by grain boundaries; the faces between two grains. Therefore, while we 
have found that premelted films at grain boundaries are 1-10 nm thick, they may dominate 
the volume of liquid through which transport is controlled. Additional relevant settings 
include how such interfacial liquidity influences the rates of atmospheric chemical processes, 
which take place on poly crystalline ice particles, the electrical conductivity of glaciers and 
ice sheets, frost heave, thunderstorm electrification and extraterrestrial ices 1 . 

In summary we have described an optical scattering study of grain boundary premelt- 
ing in ice. By doping the grain boundaries with ions to provide a colligative effect, and a 
source of repulsive screened Coulomb interactions, we found that under the experimental 
conditions such liquid enhancing effects dominate liquid suppressing long ranged attrac- 
tive dispersion forces. In all grain boundaries we find that the premelted intergranular film 
thickness increased with electrolyte concentration and at the lowest concentrations the anal- 
ysis demonstrated substantial sensitivity to the surface charge density at the grain surfaces. 
Finally, the finding of films of thicknesses up to about 10 nm has a range of immediate 
environmental applications. 
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